Rise of Kp Total Cross Section and Universality 
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Abstract 

The increase of the measured hadronic total cross sections at the highest energies 
is empirically described by squared log of center-of-mass energy yfs as otot — -Blog 2 s, 
consistent with the energy dependence of the Froissart unitarity bound. The co- 
efficient B is argued to have a universal value, but this is not proved directly 
from QCD. In the previous tests of this universality, the p(p)p, tt^p, and K^p 
forward scatterings were analyzed independently and found to be consistent with 
B pp ~ B^p ~ Bk p , although the determined value of Bk p had large uncertainty. 
In the present work, we have further analyzed forward K^p scattering to obtain 
a more exact value of Bk p - Making use of continuous moment sum rules(CMSR) 
we have fully exploited the information of low-energy scattering data to predict the 
high-energy behavior of the amplitude through duality. The estimation of Bk p is 
improved remarkably, and our result strongly supports the universality of B. 



*e-mail: mishida@wisc.edu . Address until March 2012 : Department of Physics, University of Wis- 
consin, Madison, WI 53706 



1 



1 Introduction and Summary 



The increase of total cross sections a to t in high-energy region is described [T] by squared log 
of center-of-mass energy y/s as <r tot ~ B log 2 s consistent with the energy dependence of 
Froissart unitarity bound [2]. The COMPETE Collaboration has further assumed <7 to t — 
B log 2 (s/so) to apply to all hadron total cross sections with a universal value of -B.[3] The 
universality of B was theoretically anticipated in RefsjH |5] and recently inferred from a 
color glass condensate [6] |7] of QCD. However, there is still no rigorous proof based on 
QCD. 

In the previous work[8j, to test this universality empirically, the p{p)p, 7r T p, and K T p 
forward scattering amplitudes were analyzed, and the values of B, denoted respectively as 
B pp , B np , and Bk p , were estimated independently. The resulting values were consistent 
with the universality, B pp ~ B np ~ Bk p , although there was still a rather large uncertainty 
in B Kp , B Kp = 0.354 ± 0.099mb. 

In the present work, to reduce the uncertainty of B^ p we have refined the analysis 
of forward K T p scattering. By employing continuous moment sum rules(CMSR)[HJ QUI 
[TTJ [15] for the crossing-even amplitude, we have fully exploited the information of low- 
energy scattering data to predict the high-energy behavior of cr tot through duality. We 
use two CMSRs with parameters e = —1, —3 and, following ref. [12], the estimated values 
of the cross-section integrals are treated as data points. The unphysical region of the 
integral is problematic. The pole of A(1405) in the unphysical region gives a relatively 
large contribution in K~p amplitude. We estimate its contribution by the coupled-channel 
study by A. D. Martin[13]. The cr to t an d p-ratios (=Re //Im /) of K^p forward amplitudes 
in high-energy region are fit simultaneously with two CMSR's data points. The resulting 
value is B^ p = 0.328 ± 0.045mb. The error is improved remarkably, less than half of the 
previous estimate. By comparing with previously determined values of B pp and B np , the 
universality, B Kp ~ B wp ~ B pp , is strongly supported. 



2 Formulas 

We inherit the notation and definitions from the previous workjS]: v{k) is the energy(momentum) 
of the kaon beam in laboratory system, v = \Jk 2 + m 2 K where uik is mass of K T . It is 
related to the center of mass energy \fs by 

s = M 2 + m 2 K + 2Mis . (1) 

The K T p forward amplitudes are denoted as f K ^ v {y). The total cross sections cr^t p are 
given by Im f KTp {v) = (k/An)a^ p through optical theorem. The crossing-even/odd 
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amplitudes / (±) (» are given by f {± \v) = (f K p (v) ± f K+p (u))/2. The f (+) (u) is related 
to the A'(u,t) amplitude at t = 0[T2] as / (+) 0) = A'(u,t = 0)/4vr. Real and imaginary 
parts of /^(y) are assumed to take the forms[8] 

e(+) _ JLf^^. V ■ - ^ ■ ^' ^ " 



Im/W = — C2 log 2 + Cl log + c +— (2) 

m K \ m K rriK J m K \m K J 

Im/ „ = h_(JL)"" (3) 

Re/ ,» = ^( 2c2log _^_ + Cl )_^(_i_)- cot ^ + /(+ ) ( o) (4 ) 

Re/H = ^.(JLY^L , (5) 



rriK \w,k 

in the asymptotic high-energy region, where we assume Im is given by log 2 f (and log v) 
terms in addition to the ordinary Pomeron(co) term and Reggeon(/3p/) term. Similarly 
Im /(-) is given by exchange of p, cu-meson trajectories in Regge theory with degenerate 
trajectory intercepts assumed. Their intercepts are taken to be an empirical value a P > ~ 
ay ~ 0.5 . The real parts of f&t are obtained from their imaginary parts by using 
crossing-symmetry, f^(— vr — ie) = ±f^'(u p + ie), except for a subtraction constant 
/ (+) (0). There are the total six parameters: 02,1,0, 0p',v, f^ + \0)- The C2-term in Im 
dominates cr to t in the high-energy region and it is related to Bk p by 

B Kp = 47ic 2 /m 2 K . (6) 



3 Continuous-Moment Sum Rules 

We can exploit the low-energy scattering data to predict the amplitudes in high-energy 
region by using continuous-moment sum rules(CMSR)[IT]. In ref.[12] a compact form of 
CMSR is given in ttN scattering, where Regge-pole contributions are parametrized in a 
form satisfying crossing-symmetry and convenient for CMSR. The crossing-even A'(u, t) 
amplitude was taken to be A'(v, t) = ^2 i=Ppl P ,,[— 1%{y1 — y 2 ) a ^ 2 } in the asymptotic energy 
region v > V\ . Here Vq is the normal threshold Uq = // + t/(4M), where //(M) is charged 
pion(proton) mass, and (uq — v 2 ) is analytically continued to [y 2 — v^)e~ % ' K in v > vq 
for A'(u,t) to have the correct phase factor of Reggeon-exchange amplitudes. The finite- 
energy sum rule for vA'(v,t) is given by 

r*,*uM-*)-+»»*\= E 7 ,W-g)'- + "^[fK- f -i)] (7) 

Juo i=RP'P" ai e + 1 
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where e is a continuous parameter and the nucleon pole term contribution should be added 
to the left hand side(LHS). 

If we know completely the scattering amplitudes in low-energy region from exper- 
iments, the high-energy amplitudes are predicted via CMSR through analyticity. In 
ref.[12], by using the CERN results on ttN scatterings up to 2 GeV, the low energy 
integrals on the LHS of Eq. (J7j) were evaluated in the region < — t < 1 GeV 2 and 
0(— 1) > e > — 5 for t = 0(t < 0) in steps of 0.5 . These numbers were treated as data 
points, and simultaneously fit as data to the Reggeon parametrization of the asymptotic 
high-energy region. 

This method can be applied to the analysis of forward K^p scattering. Our f^ + '(v) 
amplitude corresponds to A' /{Air) with t — 0, and the CMSR is given by 

f 1 dv vim \{m\ - z/ 2 )- (e+1)/2 / (+ V)] = [ * du vim [(m| - ^T (e+1)/2 /i s + V)]-(8) 
Jo Jo 

The LHS should be evaluated from low-energy experimental data, while the RHS is ana- 
lytically calculated by using the formulas (j2J) and (jlj). 

If we take non-odd values of e, Re f( + \v) data in low-energy region are necessary as 
inputs. However, experimental data [13] of Re f(v)(or p ratios) for K T p are poorly known 
for k < 5 GeV. In this situation we are forced to select e as odd integers; specifically we 
take e = —1, —3. 

The CMSRs with e = —1, —3 are equivalent to the n = 1, 3 moment sum rule|10j. 

- f * dv z/Tm fM(y) = - f * dv z/Tm /«(*,), (9) 

7T Jo ^ Jo 

where - is from our convention. The LHS of Eq. ([9]) is evaluated in the next section. 



4 Evaluation of integrals from experimental data 

The LHS of Eq. ([9]) is obtained by averaging the integrals of Im f K+ P and Im f K p , which 
are evaluated separately from experimental data. 

The K + p channel is exotic and it has no contribution below threshold, v < m^. Thus, 
its integral region is rrix to v\. By changing the variable from v to k, the relevant integral 
is given by 

- f * dv z/Tm f K+ ? = - I" dk kv n ^K^ (10) 

7T Jo 71 Jo 47r 
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where T\ = \Jv\ — m 2 K , which is taken to be some value in asymptotic high-energy region. 
We take V{ = 5 GeV. Actually the integral of RHS of Eq. ( FlOl) is estimated by dividing 
its region into two parts: In the low-energy part, from k — to kd, there are many data 
points[15]. They are connected by straight lines and the area of this polygonal line graph 
can be regarded as the relevant integral. In high-energy part, from kd to z/[, we use the 
phenomenological fit used in our previous work[8j. The dividing momentum kd is taken 
to be 3 GeV. As a result we obtain 



for n — 1 , 3 respectively. 

On the other hand, K~p amplitude includes Born terms of E°, A poles, which corre- 
spond to the nucleon pole term in nN amplitude. K~p is an exothermic reaction with 
open channels, An and En, below threshold which give contributions to Im / p (u) in 
unphysical regions, < v < m K and < v < m K , respectively, where v^/v^ is 
An /En threshold energy, z/ A7r / S7r = [(M A o/ 2 o + /i) 2 - M 2 - m 2 K ]/(2M). A large contribu- 
tion from 1 = A (1405) is expected to be in the unphysical region. In order to estimate 
these contributions we adopt a coupled-channel analysis by A. D. Martin[13j. In his anal- 
ysis the cross sections of K~p — > K~p, K°n, En, An ; K%p —> K®p, An + ; cr tot (K%p) , and 
Re f K P ' K n were reproduced successfully. Martin included a correction from Coulomb 
scattering following the scheme of Dalitz and Tuan[16j for channels like K~p, which gives 
a fairly large effect in very low energy region. Here we use the purely strong-interaction 
part of Martin's amplitude to estimate the unphysical region contribution. We consider 
this is the most reliable way to estimate the unphysical region. The relevant channels 
are the three 1 = 1 channels KN, En, An and the two 1 = channels KN, En. The 
S'-wave I = 1(0) scattering amplitudes 7 l/=1 (2 n/ =o~) are parametrized by the i^-matirix 
(M-matrix) as T I=1 = K(l — iqK)^ 1 (T I=0 = (M — iq)~ l ), where q is a diagonal ma- 
trix of the channel cm. momenta denoted as q = diag{p,pj],p\}(diag{p,pj]}). The real 
symmetric 3-by-3 matrix K is taken to be constant, including six parameters. The M 
matrix is taken to be effective range form, M = A + Rp 2 , which makes it possible to 
describe A(1405) resonance. Here p is the KN cm. momentum which is continued below 
threshold as p = i\p\. Real symmteric 2-by-2 matrices A, R are taken to be constant, and 
M includes six parameters. 

By using the elements of K and M, the invese of 1=1,0 KN scattering amplitudes 




{ 



20.348(41) + 72.435(184) = 692.795(189)GeV 
115.15(27) + 1294.90(3.48) = 1410.05(3.49)GeV 3 



(11) 
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Tkk'° are given explicitly [T7] by 

(tITk' )- 1 = Aj* lja +ip 

A I=1 = K KK + ^-(K K xB 2 + K KA B 3 ) 

B 1 = (1 - ip s K^)(l - ip\K AA ) + PsPaK^a 
B 2 = Pt.Pk{K K y,Kak - K KA K SA ) + Wt.K K y, 
B 3 = Pt,Pa(KkaKee ~ KksKsa) + ip\K KA 

Ajlo = M KK - (12) 
M S2 - ipy, 

where the subscript kk means (1,1) element corresponding to KN — >■ KN. Other sub- 
scripts are used similarly. Ai = ai + ibi are the S"-wave KN scattering lengths and 
M KK = A KK + R KK p 2 etc. The best-fit values of the relevant 12 parameters are given [T5] 
with errors in ref.|13j. 

Our K~p forward amplitude f K p (u) is related[19] to the T^ K by 

f^ p M = ^[T£2 + T£*]. (13) 

There are no <7 tot p -data reported below k = k s = 0.245 GeV in Particle Data Group [TB]. 
The Im f K p obtained by Martin is also utilized in this energy-region ttik < v < v s where 
v s = ^k 2 + m 2 K . 

The relevant integral is evaluated as follows: 



2 

7T 



1 z/Tm f K v {v)dv 



a 2 2 r Vs 2 r d 2 f Ul 

1§Vb,r(- M R + M + u b,r) + - dv- + - dk--- + - dk 



"An 



-0.106 ± 0.060) + (0.508 ± 0.095) + (35.481 ± 0.069) + (108.499 ± 0.405) , . 
-0.0004 ± 0.0010) + (0.110 ± 0.017) + (191.28 ± 0.54) + (1928.04 ± 7.46) ^ ' 



144.382 ± 0.426 GeV n = 1 
2119.43 ± 7.48 GeV 3 n = 3 



3 „_o (15) 



where z/ BjR = (M| — M 2 — m| r )/(2M). The 1st bracket of Eq. (JUJ) is the Born-term, 
which is estimated by using g 2 ^ = 13.7 ±1.9 and g^ = 0. The error comes from the 



6 



upper limit of g^ < 3.7 ± 1.3 [13]. The 2nd bracket, corresponding to the integral region 
v k-k < v < v s , is also estimated by Eq. ([13]) of Martin's amplitude. The error comes 
from the Rkk = 0.41 ± 0.10 fm[13], which gives the largest uncertainty among all the 
parameters. The 3rd and 4th brackets are estimated from the polygonal line graph and 
phenomenological fit [8], respectively. For n = 1 a small but sizable contribution comes 
from 1st term and 2nd term. The corresponding errors affect the final value of Eq. ([15]) 
but they are not main sources of its error. For n = 3 the 1st and 2nd terms are both 
negligible. Thus, the uncertainty from the unphysical region, especially from A(1405) 
pole, is considered to be insignificant in use of the values of Eq. ( [T5]) . 
By averaging Eqs. ([TT!) and ([15]) we obtain 



7T 







(+) . 118.589 ± 0.233 GeV 

Ulmf ^ dv = < 1764.74 ±4.13 GeV 3 (16) 



for n = 1, 3 respectively. These two values are treated as low-energy data points and they 
are fit simultaneously with the data in asymptotic high-energy region. 



5 Results and Concluding Remarks 

The a*~ p , a^ p , p K ~ p , and p K+p (more precisely]!]] Re f K ~ p , and Re f K+p ) with 
k > 5 GeV, which are given in ref. [15], are fit by using the formula in §2. The numbers of 
parameters are six: 02,1,0, $p',v-> / (0)- We fit to the CMSR data points of n = 1, 3 ( [TBI) 
simultaneously. We considered three cases: i) the case including a n = 1 datum of ( [TBI) , 
ii) the case including a n = 3 datum of ( Tl6l) . and iii) the case including both n = 1,3 
data of ([T6]) . The results are compared with iv) the case with no use of CMSR and also 
our previous analysis II[8] . 

The best fit parameters and x 2 values are given in Table [TJ 

All fits are successful ( x 2 /deg.freedom< 1 ). The best-fit \ 2 values of i),ii),iii) are 
almost the same as the case iv) with no use of CMSR, suggesting that the CMSR works 
well in the fit. Results of the best fit in the case iii) are shown in Fig. [TJand|2j 

The error of C2 in the case iii) is much smaller than the case iv), and is greatly improved 
from our previous analysis II[8] . Correspondingly the parameter B Kp associated by the 
ln 2 s dependence is given by 

Present result iii) Previous analysis II [8] 

c 2 = 0.01634 ±0.00223 <- 0.01757 ± 0.00495 (17) 

B Kp = 0.328 ± 0.045 mb <- 0.354 ± 0.099 mb 
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cr(mb) 




0.1 1 10 100 1000 GeV 



Figure 1: Results of the fit to a^ p (mb). Upper(blue) data represent K~p and 
lower(black) data are K + p. The horizontal arrow represents the energy region of the fit 
to the asymptotic amplitude. The solid curves are the asymptotic Reggeon amplitudes, 
which are extrapolated down through the resonance region. 



Table 1: Best fit parameters and \ 2 '- &tot P an d Re f KTp in k > 5 GeV are fit simulaneouly 
with i) n = 1 CMSR datum, ii) n — 3 CMSR datum, iii) n = 1 and 3 CMSR data. 
The results are compared with the case iv) with no CMSR data and with the previous 
analysis [8], where a^ p in k > 20 GeV and Re f KTp in k > 5 GeV are fit using the 
finite-energy sum rule with region of the integral 5<fc<20 GeV as a constraint. 



case 


c 2 


Cl 


c 




Pv 


/<+>(0) 


X 2 /(N D — N P — 1) 


i) n = 


1 


0.01652(224) 


-0.1241 


1.152 


0.2702 


0.5741 


1.180 


143.21/(165-6-1) 


ii) n = 


: 3 


0.01724(256) 


-0.1339 


1.188 


0.2110 


0.5736 


1.609 


143.54/(165-6-1) 


iii) n - 


= 1,3 both 


0.01634(223) 


-0.1221 


1.146 


0.2736 


0.5737 


1.178 


144.05/(166-6-1) 


iv) no 


CMSR 


0.01522(385) 


-0.1065 


1.088 


0.3726 


0.5749 


2.104 


143.04/(164-6-1) 


110 


0.01757(495) 


-0.1388 


1.207 


0.1840 


0.5684 


1.660 


63.80/(111-5-1) 



8 



p 




k 



1 10 lOOGeV 



Figure 2: Results of the fit to p K ^ p . The solid(blue) data points represent K~p and the 
open(black) data points are K + p. The horizontal arrow represents energy region of the 
fit to the asymptotic amplitude. 

There are several comments should be added: 

i) In the present analysis we take Fj~ = 5 GeV, and the <7 tot and Re / data in k > Ff are 
fit. The results are almost independent of the choice of Ff. If we take Ff = 4 GeV and 
the CMSR with e = — 1, —3 together with the data of k > 4 GeV are fit simultaneously, 
we obtain the result B^ p = 0.311 ± 0.039 mb which is almost the same result as Eq. (fT7|) . 

ii) In our previous analysis[8j, Ff = 20 GeV was taken, and the the data of a in k > 20 GeV 
and Re f(k) in k > 5 GeV were fit. The n = 1 integral with Ff = 20 GeV of LHS (1TT])(, 
which is obtained as 6802.61(10.90) GeV) is fit simultaneously to the data of the same 
energy-regions as the previous analysis [8]. The resulting value is Bk p = 0.360 ±0.110 mb, 
which is almost the same as II [8] , Br p = 0.354 ±0.099 mb so no improvement is obtained 
in this method. Hence we have adopted a different energy-region in the present analysis. 

The obtained value of B Kp together with the previous estimates of B np and B pp in 
ref. [8] are shown graphically in Fig. [3] in compared with the case with no use of sum rules. 
Our results strongly suggest the universalty of the coefficient B. It should be noted that 
the inclusion of low-energy data through the CMSR is essential in reaching this conclusion, 
especially that B^ p — B np as shown in the present work. Our conclusion is that the B 
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B(mb) 
0.6 r 



0.5 



0.4 
0.3 
0.2 
0.1 
0.0 



CMSR used 



No CMSR 



Figure 3: Values of the asymptotic parameter B(mb) for the p(p)p, vrp, and Kp total cross 
sections using CMSR(left) in comparison with the cases with no CMSR constraint (right). 
Red circles and blue squares represent the p{p)p and up, respectively, which were obtained 
in our previous analysis II [8J . Green triangles are Kp, obtained in the present analysis. 



does not depend on the flavor content of the particles scattered. 
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